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We investigate how collective modes and colored noise conspire to produce a beam halo with much 
larger amplitude than could be generated by either phenomenon separately. The collective modes 
are lowest-order radial eigenmodes calculated self-consistently for a configuration corresponding to 
a direct-current, cylindrically symmetric, warm-fluid Kapchinskij-Vladimirskij equilibrium. The 
colored noise arises from unavoidable machine errors and influences the internal space-charge force. 
Its presence quickly launches statistically rare particles to ever-growing amplitudes by continually 
kicking them back into phase with the collective-mode oscillations. The halo amplitude is essentially 
the same for purely radial orbits as for orbits that are initially purely azimuthal; orbital angular 
momentum has no statistically significant impact. Factors that do have an impact include the 
amplitudes of the collective modes and the strength and autocorrelation time of the colored noise. 
The underlying dynamics ensues because the noise breaks the Kolmogorov-Arnol'd-Moser tori that 
otherwise would confine the beam. These tori are fragile; even very weak noise will eventually break 
them, though the time scale for their disintegration depends on the noise strength. Both collective 
modes and noise are therefore centrally important to the dynamics of halo formation in real beams. 

PACS numbers: 45.10.-b, 52.25.Fi, 29.27.-a 



I. INTRODUCTION 

We recently demonstrated |]J that the combination of 
colored noise and global oscillations in intense charged- 
particle beams can create much larger halo amplitudes 
than would arise in the absence of noise. This was done 
using generic 'particle-core' models as representations 
of time-dependent potentials associated with nonequilib- 
rium beams [jy]: the 'core' established a time depen- 
dence in the form of a harmonic oscillation reminiscent of 
the presence of a global collective mode, and test particles 
orbited in response to that potential. Ever-growing halos 
were found to form despite the fact that large-amplitude 
orbits spend considerable time under the influence of the 
external focusing forces, and the frequencies associated 
with these forces differ from those associated with the 
core oscillation, a circumstance that impedes resonance. 
Thus, the noise has a key influence, boosting statistically 
rare particles to ever-growing amplitudes by continually 
kicking them back into phase with the core oscillation. 
The importance of this finding lies in the accelerator's 
extreme sensitivity to beam loss. For example, in a light- 
ion accelerator, beam impingement of just ~1 W/m at 
energies exceeding ~20 MeV will cause enough radioac- 
tivation to preclude hands-on machine maintenance |4(. 
In high-average-current machines, this amounts to just a 
few particles lost per meter, and large halos are thereby 
of practical concern, even if their outermost fringe is ex- 
tremely tenuous. 

Our previous analysis was restricted to radial orbits 
and centered on choosing the same initial conditions for 
all of the orbits. Specifically, each orbit was assigned zero 
initial velocity and the same initial radius. Because in a 
real beam each individual particle has its own distinct ini- 
tial conditions (for example, the particles would start at 
different angular coordinates), each experiences its own 



manifestation of colored noise. In other words, the noise 
was regarded to be spatially uncorrelated. Thus, we se- 
quentially computed 10,000 orbits while assigning to each 
orbit its own unique, random manifestation of the col- 
ored noise, and we cataloged the maximum amplitudes 
of these orbits. Though this approach proved sufficient 
to demonstrate the noise-enhanced production of beam 
halo, it suffers a number of shortcomings. First, it lacks 
self-consistency; with one exception, the oscillation fre- 
quencies of the core were chosen ad hoc, the exception 
relating to a space-charge-limited core. Second, because 
only a single starting radius is sampled, it lacks the statis- 
tics of a full treatment; halo particles originating from, 
e.g., different radii are excluded. Third, the contribution 
of nonradial orbits is likewise ignored. 

The present paper offers a study that, by largely cir- 
cumventing these shortcomings, is more thorough and 
systematic. Herein we consider self-consistent collective 
oscillations in the context of a general framework. Specif- 
ically, we consider a direct-current, cylindrically sym- 
metric beam and model it as a warm-fluid Kapchinskij- 
Vladimirskij (KV) equilibrium configuration. We then 
imagine the beam to be excited such that it possesses 
a self-consistent spectrum of collective, stable radial 
modes of oscillation as previously calculated by Lund 
and Davidson [£j . The associated time-dependent space- 
charge force combines with the external focusing force to 
determine the equation of motion of test particles. By 
populating the full configuration space with very many 
(typically 10 6 ) test particles, assigning each test particle 
its own random manifestation of colored noise, and then 
tracking their orbits, we compute the evolution of the 
halo. We do this for two extremes of initial particle veloc- 
ities, the first corresponding once again to purely radial 
orbits, and the second corresponding to purely circular 
orbits. The halo structure depends, of course, on (a) the 



beam parameters, which we combine into a single quan- 
tity, the space-charge tune depression, (b) the collective- 
mode parameters, specifically their amplitudes in that 
their frequencies are determined self-consistently, and (c) 
the noise parameters, specifically the noise strength and 
autocorrelation time. 

In the investigation to follow, we quantify the vari- 
ous parametric dependencies. Section [H] explains our 
methodology in detail. Section ITTT1 presents an extensive 
array of results that together quantify how the beam and 
noise parameters conspire to produce large halos. In- 
cluded is an interpretation of the underlying dynamics 
in terms of the breaking of Kolmogorov-ArnoPd-Moser 
(KAM) tori |6| due to the presence of noise. Section ITVl 
concludes by briefly summarizing the findings and, in 
view of them, identifying related phenomenology that 
will likely be inherent to fully self-consistent large N- 
body simulations of real beams. 



II. METHODOLOGY 

As our foundation, we adopt directly the formalism 
of Strasburg and Davidson Q, hereafter called SD. We 
consider an intense, direct-current charged-particle beam 
propagating in the z direction at constant speed through 
a transport channel that imposes a constant, cylindrically 
symmetric, linear transverse focusing force. The equi- 
librium beam is a warm-fluid Kapchinskij-Vladimirskij 
equilibrium, and collective modes are superposed upon 
this equilibrium. These modes correspond to stable, ax- 
isymmetric flute perturbations and derive from lineariz- 
ing the respective Vlasov-Maxwell-Poisson equations [5( • 
The influence of the beam's self-fields on particle trajec- 
tories is properly included within the framework of the 
paraxial approximation. 

We incorporate the beam parameters by way of the di- 
mensionless self-field perveance K given per the gaussian 
system of units as 
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wherein p is the line density (number of particles per 
unit length), q and m are the particle charge and mass, 
respectively, (3 and 7 are the usual relativistic factors, 
and c is the speed of light. The perveance then folds into 
the space-charge tune depression 7/ as 
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in which R is the radius of the equilibrium beam and u> / 
is the angular frequency associated with the bare external 
focusing force. This parameter lies in the range < r\ < 
1, the lower bound corresponding to the space-charge- 
limited beam, and the upper bound corresponding to zero 
space charge. 



SD tabulate the potentials and frequencies correspond- 
ing to all of the axisymmetric flute modes. These are col- 
lective normal modes, and as such are calculated using 



linear perturbation theory J5J ■ The frequency of the 
such mode is given by 
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For their studies of particle dynamics, SD concentrate on 
the two lowest-order radial modes, n— 1, 2, and we shall 
do likewise. We normalize the radial coordinate in terms 
of the radius R ; however, unlike SD, we normalize time 
t in terms of the angular frequency u>f, i.e., t—nuft. In 
effect we are setting R — 1 and ujf = 1 . 

The axisymmetric flute modes are distinctly different 
from breathing modes. The most elementary distinction 
is that the beam boundary is static (i? is constant) in 
the case of flute modes, but it oscillates in the case of 
breathing modes. In both cases the beam is root-mean- 
square (rms) mismatched; however, for the flute modes 
the beam envelope is matched whereas for the breath- 
ing modes the envelope is mismatched. We shall there- 
fore use the terms "envelope-matched" and "envelope- 
mismatched" to refer to beams with axisymmetric flute 
modes and breathing modes, respectively. For the warm- 
fluid KV beam, the equilibrium density profile exhibits 
a step-function discontinuity at the boundary. In turn, 
the flute modes likewise include a discontinuity in the 
density profile at the boundary. For example, consider 
the KV beam to be excited by the n = 1 flute mode. 
The density profile inside the beam is always uniform, 
but its magnitude oscillates. To conserve particle num- 
ber, this mode includes an oscillating surface charge, i.e., 
the density profile exhibits a Dirac delta function at the 
(stationary) envelope radius such that the integral over 
the beam volume is independent of time. By contrast 
the lowest-order breathing mode entails a self-similar os- 
cillation; the envelope radius oscillates, and the number 
density likewise oscillates but is everywhere uniform. 



A. Equation of Test-Particle Motion 

To explore the dynamics of halo formation, we com- 
pute orbits of test particles that move in the total poten- 
tial formed by the superposition of the external focusing 
potential and the space-charge potential. The test par- 
ticles contribute nothing to the total potential and do 
not interact with each other. This means we treat the 
coarse-grained form of the beam's distribution function, 
thereby ignoring, e.g., discreteness effects from the in- 
dividual point charges that comprise the beam. Using 
the formalism herein, it may be possible to mimic dis- 
creteness effects by modeling them as appropriately weak 
gaussian white noise @|, i.e., noise that has zero autocor- 
relation time, but we refrain from doing so in favor of 
concentrating on the influence of colored noise. 

The equation of test-particle motion decomposes into 
two regimes, one for which the normalized radial coor- 



dinate r < 1, and the other for which r > 1. If only the 
n = 1, 2 normal modes are excited, then the SD equation 
of test-particle motion with our normalization is 



+ r\ r — 



L 2 



- (i - ??>{vri cosM?r )t] 



r 2 (l - -r 2 )cos[w 2 (?7 2 )t]} = for r<l; 



L 2 



1-rj 2 



= 



forr>l; (4) 



in which L is the dimensionless angular momentum, T n is 
the ratio of the rms electrostatic energy contained in col- 
lective mode n to that contained in the equilibrium beam, 
and ui n is given by Eq. © after setting ujf = 1 . The con- 
stants L, Ti, and r 2 may be regarded as free parameters, 
the former in regard to the 'geometry' of the test-particle 
orbit, and the latter in regard to the 'amplitudes' of the 
respective collective modes. However, because it derives 
from linear perturbation theory, for Eq. (0J to be valid, 
both Ti and r 2 must be small compared to unity. Note 
that the tune depression r\ manifests itself not only in the 
mode frequency, but also in the frequency characterizing 
the effective focusing force acting on the test particle. 
Hence, any noise that shows up in the tune depression 
influences both of these frequencies. 



B. Colored Noise 

It is at this point that we depart in an important way 
from the SD treatment, for we wish to assess the extent to 
which noise, in combination with the collective mode(s), 
influences the particle dynamics. This is a problem of 
practical importance; noise is unavoidable in real acceler- 
ators because they are imperfect. Machine errors, as well 
as transitions, will feed space-charge fluctuations in that 
the beam evolves self-consistently in response to external 
influences. Examples include forces from image charges 
due to irregularities in the accelerator hardware as well 
as radiofrequency and magnetic field errors, and in the 
lab frame the errors may themselves be time-independent 
or fluctuating due to, e.g., jitter |9(. ^From the perspec- 
tive of a beam particle, the effect of all of these machine 
imperfections is to impart time-dependent noise on the 
particle orbit, and thus we seek now to include this noise 
in the equation of test-particle motion. Concerning our 
upcoming analysis, for zero noise we of course reproduce 
the dynamics that SD describe. Thus, any differences 
that show up with nonzero noise are attributable solely 
to the presence of the noise itself. Our main interest is 
to quantify how this noise influences the process of halo 
formation, and do so to an extent well beyond what we 
did previously. 

Following the philosophy and procedure of our earlier 
investigation |lj , we add gaussian colored noise that sam- 
ples an Ornstein-Uhlenbeck process [l(| ■ We do so in 
terms of a frequency fluctuation Suj(t). Because the tune 
depression rj incorporates the space charge, we define this 



frequency fluctuation in terms of a fluctuating tune de- 
pression in a manner consistent with Eq. JSJ: 
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The frequency fluctuation 5oj{t) henceforth represents 
the noise. Thus, everywhere it occurs in the equation of 
motion, Eq. Q, the quantity rj 2 is replaced by r] 2 +Si] 2 (t), 
with Srj 2 (t) given by the last expression in Eq. JSJ) above. 
Note, for example, that the noise will still manifest itself 
in Eq. (0} even if no collective mode is excited. 

The first two moments of 5ui(t) fully determine the 
statistical properties of the noise: 
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in which t c denotes the autocorrelation time, i.e., the 
time scale over which the signal changes appreciably. The 
special case of white noise corresponds to the limit t c ~^0. 
After generating a colored-noise signal using an algorithm 
first presented in Ref. [ll|, we compute \A\ <-> (I'M) 
which then constitutes the measure of noise strength. 
Example manifestations of colored noise for various noise 
strengths and autocorrelation times are plotted in Fig. 1 
of Ref. 0. 

The noise should be viewed from the perspective of 
the charged particle as it progresses along its trajectory. 
A typical particle will respond to space-charge fluctua- 
tions that change over a time scale comparable to, e.g., a 
plasma period (which is thus a measure of the minimum 
autocorrelation time). It will likewise respond to stochas- 
tic changes in the external fields arising as the particle 
transits the beam line. Such changes can correspond to, 
e.g., the spacing between hardware components. Accord- 
ingly a hierarchy of autocorrelation times comprises the 
actual noise a particle experiences. In a real beam each 
individual particle will have its own distinct initial con- 
ditions and thus experience a manifestation of the noise 
differing from that seen by each of the other particles. 
Hence, at each successive time step during an orbit in- 
tegration, a randomly generated frequency fluctuation is 
computed in keeping with the specified statistical proper- 
ties of the colored noise. Then this frequency fluctuation 
is converted to a fluctuation in the space-charge tune de- 
pression per Eq. JSJ) , and the so-modified tune depression 
is inserted into the equation of motion, Eq. (0J, whereby 
the fluctuation influences the next time step. How the 
noise quantitatively affects halo formation depends on 
its strength and its autocorrelation time, dependencies 
that we quantify herein. 



C. Initial Distribution of Test-Particle Orbits 

In keeping with the objective of retaining as much real- 
ism as possible, we choose a beam intermediate between 



the space-charge-limited beam (77 = 0) and the zero-space- 
charge beam (rj = l). Specifically, we keep r\ fixed at the 
value 1] — 0.3 for our entire investigation. 

We choose an initial distribution of test particles that 
spans all of the dynamically interesting regions of con- 
figuration space and falls gradually to a low-density tail. 
Of course there are numerous ways to do this; one is 
to choose a distribution corres ponding to a configuration 
of thermal equilibrium (TE) J12, ll3J |- We construct a 
cylindrically symmetric TE configuration of test charges 
following a procedure recentl y us ed to devise spherically 
symmetric TE configurations 14] . The associated dimen- 
sionless Poisson equation is 
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wherein R is a dimensionless radial coordinate, n(R) 
denotes the number density normalized to the central 
density, fi is a dimensionless quantity governing the 
strength of the external focusing force vis-a-vis the col- 
lective space-charge force, and $(i?) is the dimensionless 
space-charge potential. For the value of VL we choose 
il= (l + l(T 3 - 5 -l(r 9 )/\/2. Integrating Eq. numer- 
ically using this specific value of fi yields a TE density 
profile corresponding to a tune depression ?/~0.3 13]. 

Length and time are normalized differently in Eq. ({Y| 
than in Eq. Q). In keeping with the desire to span all 
of the dynamically interesting regions of configuration 
space, we simply rescale the density distribution n(R) 
calculated from Eq. J7|) so that its rms radius R corre- 
sponds to the full radius R of the warm-fluid KV dis- 
tribution: R = R = 1. This clearly places a sizeable 
population of test particles, that corresponding to much 
of the density tail, outside the KV 'core'. It also mim- 
ics, e.g., an inference from the recent beam-halo exper- 
iment at the Low-Energy Demonstration Accelerator at 
Los Alamos National Laboratory that the input beam 
for this experiment carried a sizeable tail in its distribu- 
tion [TJI. For all of our investigations the initial radii 
of the test particles follows this distribution. Most of 
our simulations involve N = 10 6 test particles, a num- 
ber sufficient to constitute a good statistical sample. In 
principle, the tenuous tail of the density profile extends 
to infinity, but in practice there is a finite radius to the 
A-body representation of the density because N is finite. 

For most of our investigations the initial test-particle 
velocities are all set to zero, corresponding to purely ra- 
dial orbits, in which case we then replace r(t) by x(t), 
and r(t) < 1 or > 1 by \x(t)\ < 1 or > 1, respectively, in 
Eq. J2J. We also, however, consider another 'limiting' 
case, that for which all the orbits are initially circular. 
Given a radius r c of the initially circular orbit, the re- 
spective dimensionless angular momentum L, a quantity 
taken to be conserved, follows from Eq. (0J: 
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As is shown and discussed in Sec. IIII C 31 below, the in- 
fluence of noise on circular orbits that start with r c < 1 is 
essentially the same as for the purely radial orbits. Con- 
sequently, the halo population is similar for both cases. 



D. Orbit Integrations 

We integrate the equation of motion using a fifth-order 
Runge-Kutta algorithm with variable time step [lq tak- 
ing the initial time step to be 0.01 'differential-equation' 
(DE) units. We evolve each orbit for a total time 512 
DE units, which corresponds to 40-60 orbital periods de- 
pending on the initial conditions for the respective orbit. 
Thus, for example, the total integration time is compara- 
ble to the transit time of the beam through a large pro- 
ton linear accelerator such as that associated with the 
Spallation Neutron Source [Tjfl • In the absence of time- 
dependence and noise, i.e., with Ti = T 2 = (\6v\) = 0, 
the algorithm conserves energy within a fractional error 
10~ 9 at each time step and within 10~ 7 over the whole 
integration. 

Our investigation spans a broad sector of the param- 
eter space in that the space-charge tune depression, set 
at i] = 0.3, is the only parameter that is never varied. We 
treat all combinations of the following parametric val- 
ues: mode amplitudes r^ =0.05, 0.10, and 0.20; noise 
strengths (\Suj\} = 0, 0.001, 0.01, and 0.1; autocorrelation 
times t c = 0.5, 1, 2, 10, 80, and 160; and test-particle 
sample sizes N=10 2 , 10 3 , 10 4 , 10 5 , and 10 6 . Most of the 
plots shown herein pertain to the specific choice t c = 80; 
however, excursions to lower and higher values are in- 
cluded to provide a check on the sensitivity of halo for- 
mation to the autocorrelation time of the noise. In a real 
machine, of course, the noise will incorporate a range, or 
'superposition', of autocorrelation times and strengths. 
In addition, with one exception (Fig. [SJ), all of the plots 
we show pertain to radial orbits, i.e., cases for which all 
of the test particles have L = in Eq. (gj . 



III. INFLUENCE OF COLORED NOISE ON 
HALO FORMATION 

SD explored the dynamics of test-particle motion in the 
absence of noise. The form of their equation of motion 
differed slightly from our Eq. Q in that they normalized 
the time in terms of the space-charge-depressed focusing 
frequency rather than the external focusing frequency. 
Notwithstanding the different normalization, the physi- 
cal content remains unchanged. SD discovered that the 
time-dependent potential associated with the presence of 
a collective axisymmetric flute mode, even if only weakly 
excited, establishes a chaotic region of phase space in the 
outer regions of the beam. They also found that this 
feature is not present in the phase space of an envelope- 
mismatched beam having a similar level of rms mismatch 
but no collective mode (cf. Fig. 9 of Ref. 0|). An impor- 



tant consequence of the chaotic sea is that orbits enter- 
ing it can stochastically explore a larger region of phase 
space, thereby gaining more energy and correspondingly 
larger orbital amplitude. SD consequently demonstrated 
that the excitation of collective modes leads to a halo 
significantly larger than that generated by an envelope 
mismatch, and that the difference is due to destabiliza- 
tion of KAM surfaces by the collective mode. As we will 
now show, the presence of colored noise substantially en- 
hances the influence of the collective modes. 



A. Orbital Dynamics 

Halo formation is inextricably linked to the dynam- 
ics of individual orbits. Consequently, a close inspection 
of what happens to an individual orbit because of the 
noise will be instructive. We arbitrarily select an orbit 
that originates deep in the interior of the beam; the ini- 
tial conditions are x = —0.733407 and x = 0. We then 
integrate the orbit for 2048 DE units to obtain good fre- 
quency resolution in its power (Fourier) spectrum. The 
trajectory and power spectrum of this orbit are plotted 
in Fig. [2 with Tj = 0.1, T 2 = 0, and t c = 80, and for a 
sequence of successively increasing noise strengths. With 
zero noise the power spectrum is sharply peaked at a sin- 
gle frequency indicating that the orbit is periodic, hence 
regular, and the trajectory x{t) clearly reflects this pe- 
riodicity. However, in the presence of even weak noise 
the orbit clearly becomes chaotic, having a power spec- 
trum that features continua. A useful measure of chaos 
is the number of frequencies Kf that together contain a 
given fraction / of the spectral power. This measure is 
called the 'complexity' of the orbit; a common choice is 
/ = 0.9 [l4|. Accordingly, a broader spectrum indicates 
a higher degree of chaos. Concerning the orbit in Fig.^ 
the noise has obviously placed it in a chaotic region of 
the phase space established as a consequence of the time- 
dependent potential associated with the collective mode. 
How this happens is clarified in Sec. IIII El below. 

As is also apparent from Fig. ^ the degree of orbital 
chaoticity as quantified via the complexity Kqq is not 
necessarily a simple, i.e., monotonic, function of the noise 
strength. To reiterate, the power spectrum, hence the 
complexity, derives from the history of the orbit and 
thereby reflects a superposition of successive short-time 
behaviors. An orbit that spends a relatively large fraction 
of time at large amplitudes, over which the net force is 
predominantly that of the harmonic external potential, 
will tend to be 'more regular' and have smaller Kq.q. 
As concerns a single, specific orbit (so no phase-space 
statistics are involved), what matters is not so much the 
amplitude of the noise, but rather whether a sequence 
of noise-induced kicks happens to make the orbit more 
chaotic, and these kicks are, of course, unpredictable a 
priori. Our experiments indicate that a sequence of kicks 
leading to increased orbital chaoticity and/or increased 
orbital amplitude will occur sooner for some orbits and 



later for others. In simulations involving many test parti- 
cles distributed over a range of initial conditions, features 
of the evolving test-particle distribution are thus mani- 
festly statistical. 



B. Evolution of the Halo Amplitude 

We now evolve initial distributions of N = 10 6 test 
particles constructed per the prescription of Sec. Ill Cl As 
the orbit integrations progress, we record a 'snapshot' of 
the test-particle positions once every eight DE time units. 
This interval approximately corresponds to the period of 
a typical orbit in the unperturbed SD potential, which we 
call the 'dynamical time' to'- 8 DE units ~1 tr>. For every 
snapshot we record the largest radius reached by any of 
the N particles; the collection of these radii represents 
the evolving halo amplitude Ru{t)- 

Example results are plotted versus time in Fig. [21 for 
which the mode amplitudes and autocorrelation time are 
fixed at Ti = 0.05 or 0.1, T 2 = 0, and t c = 80, and the 
noise strength (|<W|) is varied from zero upward. As the 
figure indicates, in the absence of noise the halo ampli- 
tude is quasiperiodic, and its time-averaged value stays 
the same, i.e., it does not grow. This is as expected 3, 7]. 
A particle 'resonantly' coupled to the collective mode is 
kicked to larger amplitudes. However, because its or- 
bital frequency changes as its amplitude changes, at suf- 
ficiently large amplitude the particle decouples from the 
mode and its amplitude ceases to grow. Differences be- 
tween the external focusing force and the collective space- 
charge force thus impose a hard upper bound on the halo 
amplitude. The presence of noise, however, drastically 
changes this scenario. Occasionally successive kicks from 
the noise will happen to be 'just right' to keep a particle 
in phase with the mode for an effectively longer time and 
thereby push it beyond the upper bound (i.e., outer KAM 
torus) characterizing the noise-free case. The halo am- 
plitude Rh (t) continues to grow, and the growth appears 
to be almost linear with time (at least after the first few 
oscillations). Over the range of noise strengths (\Su\) we 
explore, both stronger noise and larger mode amplitudes 
enhance halo growth. Moreover, when the noise is strong 
(e.g., (\Slu\) =0.1, i.e., roughly 10% of the collective-mode 
frequency), pronounced halo growth occurs in just a few 
(~ 5) dynamical times. 

Consider the largest orbital amplitude reached by any 
particle over the course of a simulation, i.e., the largest 
halo amplitude, and denote this amplitude as max(RH)- 
This quantity will of course vary with the number of test 
particles N in the simulation. Because the number of 
particles that can be incorporated into iV-body simu- 
lations is inherently limited by available computational 
power, it is of interest to know how sensitive the halo 
amplitude can be to the choice of N. To quantify this 
sensitivity we adjust the test-particle population between 
10 2 < N < 10 6 particles. Then we perform a number of 
experiments with different noise strengths and different 



values of Ti (with T2 = and t c = 80). Results pertaining 
to r!=0.05 and 0.1 appear in Fig. Because there is a 
hard upper bound to the halo amplitude in the absence 
of noise, max(RH) is essentially independent of TV for the 
case (\6oj\) = provided TV is sensibly large. This is not 
true when noise is present; the noise establishes a quasi- 
logarithmic dependence oimax(RH) on TV, a finding that 
is in keeping with our earlier results |l|. Larger values 
of Ti yield larger values of max(RH), but the scaling of 
max{R}i) with TV remains roughly the same. 

Results presented thus far correspond to a single au- 
tocorrelation time t c — 80 ~ 10 to- What happens if t c 
is much shorter or much longer? Plots of maximum halo 
amplitude max(RH) versus t c for a sequence of increasing 
noise strengths (|<5w|), and with only the n=l collective 
mode active, appear in Fig.0J Data points in these plots 
each correspond to a sample of TV = 10 test particles; 
simulations with TV= 10 6 are found to give similar results 
but, of course, they involve much longer run times. In 
most cases the halo extent is seen to be only weakly de- 
pendent on autocorrelation time. The exception pertains 
to large noise strength: (\Sui\) =0.1, i.e., a 10% fluctua- 
tion in the collective-mode frequency, generates substan- 
tially larger halo for t c < 100. The presence of such large 
noise would seem to be anomalous in a real accelerator, 
and one might thus presume the halo amplitude will nor- 
mally be independent of hardware details associated with 
the establishment of noise correlations. However, under 
circumstances that lead to a turbulent beam as might 
reside, for example, at and just downstream of the beam 
source and at large hardware transitions, one might in- 
deed expect the particle orbits to experience large noise 
from space charge locked in the turbulent eddies. Such 
circumstances would seem normally to be transient, with 
the large-scale turbulence mixing away in a few dynam- 
ical times. Nevertheless, because it would form rapidly, 
a sizeable halo would likely arise as evidenced from the 
(|£u;|)=0.1 curves in Fig.0 



that measurably increases the halo extent beyond that 
corresponding to the n = 1 mode acting alone. As seen 
from Fig. El adding the n = 2 mode does modify the dis- 
tribution, though its effect appears to be modest. 



2. Mixing and Halo Formation 

To visualize how orbits mix, we integrate collections of 
1600 test-particle initial conditions clumped into tightly 
localized regions of phase space. The integrations are 
done both without and with noise; the evolution of these 
collections is depicted in Fig. Rows (a)-(c) pertain 
to the absence of any collective mode, i.e., Ti = T2 = 0. 
In the absence of noise [row (a)] all mixing is due to a 
frequency spread across the initial clumps arising from 
the nonlinear net force. Separation of the initially local- 
ized particle trajectories then proceeds as a power law 
in time; this is regular phase mixing, i.e., linear Landau 
damping. Noise [rows (b) and (c)] influences the effec- 
tive focusing force acting on the test particles, but this 
influence generates no significant spreading to large or- 
bital amplitudes. Turning on the n — 1 collective mode 
changes the situation completely. With Ti = 0.05, but 
in the absence of noise [row (d)] , the clumps still spread 
only to a restricted region of phase space; however, the 
same is clearly not true when noise is included [rows (e) 
and (f)]. Noise causes a far more efficient mixing. Even 
moderately weak noise, e.g., (\8u>\) =0.01, thoroughly and 
exponentially mixes particles, regardless of their starting 
points, into all regions of the phase space accessed by the 
beam. Their exponential separation into global regions 
of phase space is the principal signature of chaotic mix- 
ing [l4[. Increasing Ti further accentuates this chaotic 
mixing and causes the orbits to fill an even larger phase- 
space area. 



3. Circular vs. Radial Orbits 



C. Evolution of the Test-Particle Distribution 

1. Halo Density 

Not just the amplitude of the halo is of interest, but 
also its density profile. A convenient and meaningful rep- 
resentation of the test-particle distribution is obtained 
by calculating the percentage of particles lying outside 
a radius R; let us call this P(r > R) while noting 
P{r>R)^ 100% as i?->0. Plots of log 10 [P(r > R)] ver- 
sus R computed at £ = 512 appear in Fig. |SJ Here again, 
the mode amplitudes and autocorrelation time are fixed 
at Ti = 0.05 or 0.1, T2 = 0, and t c = 80, and the noise 
strength (|<5o;|) is varied over a considerable range. As a 
general trend the distribution spreads to larger radii, i.e., 
the halo amplitude grows, as the noise strength increases. 

One might anticipate that the n = 2 mode would cou- 
ple to a statistically small set of particles in a manner 



Thus far all simulations have pertained to radial test- 
particle orbits. Might the results be substantially differ- 
ent for orbits with nonzero angular momentum? To an- 
swer this question, we now consider the other extreme, 
that for which all test particles are on circular trajecto- 
ries at t = 0; the corresponding values of angular momen- 
tum L are assigned according to Eq. (JSJ). We then solve 
Eq. (0} without and with noise in the space-charge tune 
depression. This, however, means that we refrain from 
adding noise to the azimuthal motion. To do otherwise 
would vitiate using angular momentum as an integral of 
the motion and thereby lengthen the computations, all 
for the 'benefit' of incorporating no fundamentally new 
or important additional phenomenology. 

In Fig. [S] results for the halo amplitude Rh (t) [panel 
(a)] and halo distribution P(r > R) [panel (b)] are jux- 
taposed against those pertaining to purely radial orbits. 
Although the curves are not identical, neither are they 



systematically different. We attribute the differences to 
statistical fluctuations caused by the random noise in- 
cluded in the simulations. This finding is interesting in 
that particles on circular orbits with r c (0) > 1 essentially 
lie outside the influence of the time-dependent poten- 
tial arising from the collective modes, whereas particles 
on radial orbits do not. That halo profiles correspond- 
ing to radial versus azimuthal orbits are similar therefore 
suggests that particles initially in the 'core', for which 
r(0) < f , are the ones that dominate the process of halo 
formation independent of their initial conditions in ve- 
locity space. By reproducing Ru{t) with the initial radii 
truncated at r(0) = r c (0) = f , we verified that this sug- 
gestion is indeed true. 



D. Collective Modes vs. Envelope Mismatch 

The preceding results have illustrated how collective 
modes, orbital chaoticity, and noise in an envelope- 
matched beam collaborate to drive an ever-growing 
halo. Consider, by contrast, a beam that is envelope- 
mismatched, i.e., one that exhibits the lowest-order 

breathing mode [18|, [H| • The corresponding equation of 
test-particle motion is similar to that governing a beam 
with a single excited collective mode, except now the 
'core radius' R = R(t), i.e., the radius defining the 'in- 
side' of the beam, is a function of time. Specifically, the 
dimensionless equation governing the core radius is 



n 2 1 ~ n 2 

R+R -h--tr = ° 



(9) 



which then folds into the dimensionless single-particle 
equation of motion 



f - ri 2 
/• + /• T77r—r = for r<R(t) 



R 2 

1-t] 2 



= for r > R(t) 



(10) 



We define the envelope-mismatch parameter M = 
R(t — 0)/R o , with R = \ denoting the radius of 
the matched beam. Results presented for envelope- 
mismatched beams derive from integrating orbits in the 
coupled Eqs. @ and JTUJ. 

To obtain a rough correspondence between the ampli- 
tude Ti of the n = 1 collective mode in the envelope- 
matched beam and the envelope-mismatch parameter 
M, we imagine M to be small. We then set R —* 
l + (M-l)cos[wi(ry 2 )i], and in turn put R~ 2 (t) -*■ 1- 
2(M — 1) cos[ix>i(r] 2 )t} in Eq. (|TU|l . By comparing the end 
result with Eq. (0J (with L set to zero) and otherwise 
neglecting the different definitions of core radii, we then 
infer ri<->4(M-l) 2 . 

Now, taking the core radius to be r = R(i) as pertains 
to an envelope mismatch rather than r = 1 as pertains 
to collective modes in an envelope-matched beam has a 
profound effect on the mixing, hence the halo dynamics. 



Figure [7| [row (g)] illustrates the mixing in an envelope- 
mismatched beam. Here the mismatch is M = 1.1118, 
for which the counterpart collective-mode amplitude is 
Ti = 0.05. Accordingly, the parameters of Fig. |7Jg) are 
analogous to those of Fig. Cff). What is striking is how 
comparatively constrained the phase mixing and halo 
growth turn out to be in the envelope-mismatched beam. 
The underlying dynamics are clearly different. How so? 
The answer lies in the Poincare surfaces of section 
(PSS); as SD point out pj, these are distinctly different 
for the two cases. The PSS for the envelope-mismatched 
beam exhibits robust, densely packed KAM tori in the 
region of phase space exterior to the beam. This is true 
even if the mismatch is large. These tori inhibit the par- 
ticles from gaining significant energy and reaching large 
amplitudes. The n = 1 collective mode in the envelope- 
matched beam, by contrast, weakens the tori in the vicin- 
ity of the beam edge r = 1 . As a consequence, particles are 
then freer to move; they can stochastically and rapidly 
explore a large region of phase space. SD point out that 
the orbital amplitudes of those particles can rapidly in- 
crease as the amplitude of the collective mode is raised, 
whereas for the envelope-mismatched beam, test particles 
gain negligible energy as the mismatch is raised. SD's 
findings pertain to zero noise; we find nonzero noise sub- 
stantially magnifies them by further weakening and/or 
breaking the tori. When the KAM tori are broken, a se- 
ries of small, successive kicks can much more easily push 
a particle to ever increasing radii. Moreover, the asso- 
ciated time scale is short; significant extended halo can 
form in just a few dynamical times, i.e., orbital periods. 



E. Noise-Induced Breakdown of Tori 

To visualize noise-induced disintegration of tori with 
consequent halo formation, we plot the Poincare sections 
of 18 test-particle orbits having initial conditions that 
collectively represent the whole of configuration space. 
We take only the n = 1 mode to be excited, with am- 
plitude Ti =0.05, and integrate the 18 trajectories for a 
total time £ = 2048 (about 250 to). We do a series of these 
experiments, starting with a noise strength (|(5u;|) = 10~ 6 
and successively increasing it to 10~ 5 , 10~ 4 , 5 x 10~ 4 , 
and 10~ 3 ; for every experiment we set £ c = 80. We record 
the positions and velocities of the particles at every pe- 
riod T=2tt/u>i. For each experiment the respective PSS 
is shown in Fig.[5J different colors denote different orbits. 

The Poincare sections clarify the underlying micro- 
scopic dynamics. As the noise strength is raised, 'inter- 
nal' (lower-energy) tori are clearly the first, and thus the 
easiest, to break. With stronger noise the outermost tori 
break as well. Note that the strongest noise considered 
is only a 0.1% fluctuation of the mode frequency, and yet 
this noise breaks all of the tori [panel (f )] . It is important 
to remember that these plots lack statistical significance 
since only 18 orbits are represented. In a statistically im- 
portant, i.e., much larger, sample some number of parti- 
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cles may conceivably break through the outer tori even 
with very small noise. What the plots suggest is that this 
number should increase as the noise strength increases, 
in keeping with what one would expect intuitively. The 
'disintegration time scale' (delineating the onset of halo 
formation) is as indicated in the plots of Rh discussed 
earlier, and these plots were developed with meaningful 
statistics, i.e., with 10 6 test particles. 

We now repeat the same investigation for the cor- 
responding envelope- mismatched beam, i.e., M = 1 + 

1/2 

Iy /2= 1.1118. Here, we record the positions and veloc- 
ities of the particles at every period at which the core has 
its minimum radius; the results appear in Fig. 1101 Panel 
(a) depicts the PSS with zero noise, whereas panel (b) 
depicts the PSS with noise having the same parameters 
((|<5u;|) = l(r 3 and t c = 80) as in Fig.GJf). Although with 
this noise the KAM tori become slightly fuzzier, they 
are not yet broken, and this stands in stark contrast to 
the situation wherein a collective mode of similar am- 
plitude is active in an envelope-matched beam. Even 
when the noise has very large amplitude [(|<W|) = 10 _1 
in Fig. HOf c)]. the beam boundary in phase space ap- 
pears still to be sharply defined, although the tori have 
now obviously broken to a certain degree, particularly 
in the beam's interior. Hence, although the tori of the 
envelope-mismatched beam are seen to be robust, they 
are not indestructible; sufficiently strong noise will even- 
tually break them. 



IV. DISCUSSION AND CONCLUSIONS 

The foregoing results illustrate that collective modes 
can have a critical impact on halo formation by destabi- 
lizing the phase space near the beam boundary. Particles 
then are freer to roam by interacting with the collective 
modes, extracting energy from them, and thereby pop- 
ulating a halo. By keeping a statistically small number 
of particles in phase with collective modes, colored noise 
contributes toward not only populating the halo, but also 
expanding its extent, and it does so rapidly. 

One might reasonably question, because they are un- 
realistic, whether density discontinuities inherent to col- 
lective modes in the warm-fluid KV beam might vitiate 
our findings by imposing correspondingly unrealistic dy- 
namics. The answer would seem to be no because we 
showed in earlier work (Ref. [Jj) that adding noise to 
a perturbed thermal-equilibrium beam, a beam that is 
devoid of discontinuities and wherein the perturbation 
mimics the presence of a global collective mode, yields or- 
bital amplitudes entirely consistent with those computed 
herein. 

These matters are of practical importance to the evolu- 
tion of real beams. Transitions in an accelerator will give 
rise to various mismatches that move the beam away from 
equilibrium. Subsequent charge redistribution will trig- 
ger a hierarchy of collective modes. Unavoidable irregu- 
larities in the beamline will impose a spectrum of colored 



noise that adds self-consistently to the time-dependent 
potential associated with the collective modes. Conse- 
quently, the phenomenology that we uncovered will arise, 
as will the consequential growth of the beam's phase 
space in general, and beam halo in particular. Account- 
ing for these details therefore becomes imperative, par- 
ticularly in regard to designing accelerators for the pro- 
duction of high-average-current beams. 

Although by working with the warm-fluid model of 
a beam we have endeavored toward a treatment that 
is realistic, yet still generic, our treatment neverthe- 
less retains some shortcomings that need to be recti- 
fied in future work. These include the following: (1) 
The distribution of collective modes will evolve in a real 
beam 20]; modes will tend to dissipate in conjunction 
with the redistribution of the free energy they contain, 
a dynamic that we have neglected. However, the time 
scales over which large-scale collective modes dissipate 
are not yet well quantified, and evidence from numer- 
ical simulations suggest they may persist for hundreds 
of dynamical times |2l|. To the extent this proves true, 
our analysis reveals the attendant impact on halo forma- 
tion. (2) A real beam contains no test particles; all of 
the particles interact with one another. (3) A real ac- 
celerator will present a spectrum of colored noise, i.e., a 
distribution of noise parameters, in keeping with the ac- 
tual hardware and field irregularities. The totality of this 
phenomenology can be incorporated only by way of care- 
ful self-consistent A-body simulations that reflect both 
accurate boundary conditions and statistically accurate 
initial conditions, as well as faithfully reproduce the hi- 
erarchies of spatial and temporal scales intrinsic to the 
evolving beam. 

We have endeavored to show clearly and convinc- 
ingly that details can be important to the evolution of 
a charged-particle beam under the influence of space 
charge, in that they can make a substantial impact on 
the macroscopic evolution of its phase-space distribution. 
Accordingly, these details merit careful study. A seem- 
ingly probable outcome would be that the proper way 
to picture gcncrically a noncquilibrium beam subject to 
self-forces is in terms of an increasingly well-mixed and 
continually growing phase space as opposed to a phase 
space in which tori largely partition, and hence constrain, 
the motion of the constituent particles. This is espe- 
cially true considering that the results herein pertain to 
1.5-dimensional beams (the half dimension correspond- 
ing to time), whereas real beams are higher-dimensional 
systems, and thus their phase spaces are inherently less 
hospitable to barriers in the form of tori and cantori. 

In a paper that appeared subsequent to the submis- 
sion of our manuscript, Gerigk shows that modest statis- 
tical errors in the focusing gradient generate continually 
growing halo [22j. His initial study centers on a cylin- 
drical beam that, in the absence of these errors, is in 
equilibrium, e.g., it is both rms-matched and envelope- 
matched. He finds that the focusing errors excite oscilla- 
tions in this matched beam that then transfer energy to 



single particles. He also demonstrates that the same phe- 
nomenology applies in a more realistic three-dimensional, 
i.e., bunched, beam. Inasmuch as Gerigk was unaware of 
our earlier paper (Ref. [lj) [23, and we were likewise 
unaware of his related activity, through complementary 
investigations we have all independently arrived at the 
same conclusion: noise constitutes a continual source of 
enhanced halo production. 



Acknowledgments 



This work was supported by the Department of Educa- 
tion under Grant No. P116Z010035, and by the Depart- 
ment of Energy under Grant No. DE-FG02-04ER41323. 



[1] C. L. Bohn and I. V. Sideris, Phys. Rev. Lett. 91, 264801 

(2003). 
[2] Y.-J. Chen et al., in Proceedings of the 1991 Particle Ac- 
celerator Conference, edited by L. Lazema and J. Chew 
(IEEE, Piscataway, NJ, 1991), p. 3100. 
[3] R. L. Gluckstern, Phys. Rev. Lett. 73, 1247 (1994). 
[4] R. Jameson, Fus. Eng. Des. 32-33, 149 (1996). 
[5] S. Lund and R. C. Davidson, Phys. Plasmas 5, 3028 

(1998). 
[6] A. J. Lichtenberg and M. A. Lieberman, Regular and 

Stochastic Motion (Springer- Verlag, Berlin, 1983). 
[7] S. Strasburg and R. C. Davidson, Phys. Rev. E 61, 5753 

(2000). 
[8] H. E. Kandrup, I. V. Sideris, and C. L. Bohn, Phys. Rev. 

ST Accel. Beams 7, 014202 (2003). 
[9] J. Qiang et al., Nucl. Instrum. Methods Phys. Res. Sect. 
A 457, 1 (2001). 
[10] N. G. van Kampen, Stochastic Processes in Physics and 

Chemistry (North Holland, Amsterdam, 1981). 
[11] I.V. Pogorelov and H.E. Kandrup, Phys. Rev. E 60, 1567 

(1999). 
[12] M. Reiser, Theory and Design of Charged Particle Beams 

(Wiley, NY, 1994), cf. Section 5.4. 
[13] N. Brown and M. Reiser, Phys. Plasmas 2, 965 (1995). 
[14] C. L. Bohn and I. V. Sideris, Phys. Rev. ST Accel. Beams 

6, 034203 (2003). 
[15] J. Qiang, et al., Phys. Rev. ST Accel. Beams 5, 124201 

(2002). 
[16] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes in C (Cambridge University 
Press, Cambridge, 1995). 
[17] Spallation Neutron Source Report No. 100000000- 

PL0001-R10 (unpublished). 
[18] J. S. O'Connell, T. P. Wangler, R. S. Mills, and K. R. 
Crandell, in Proceedings of the 1993 Particle Accelerator 
Conference, edited by ST. Corneliussen (IEEE, Piscat- 
away, NJ, 1993), p. 3657. 
[19] T. P. Wangler, K. R. Crandall, R. Ryne, and T. S. Wang, 

Phys. Rev. ST Accel. Beams 1, 084201 (1998). 
[20] C. L. Bohn, Phys. Rev. Lett. 70, 932 (1993). 
[21] I. Haber, et al., Nucl. Instrum. Methods A519, 396 

(2004). 
[22] F. Gerigk, Phys. Rev. ST Accel. Beams 7, 064202 (2004). 
[23] F. Gerigk, private communication. 



10 



4.0 



x 0.0 



-4.0 



a 



iill'll'I'll'illlfl'illlllllillWlf 

!S||ijl|i|llll|il||l|l|IJI|!|i|||llllll!Jl|il 




ao VHHH 







500 1000 
time 



4.0 



0.0 - 



-4.0 
4.0 



0.0 



-4.0 
4.0 



-4.0 
4.0 

0.0 - 
-4.0 





-4-2 2 



FIG. 1: Plots of an example orbit having initial conditions 
x(0) — —0.733407, x — 0, in the presence of various noise 
strengths with I\ = 0.1, T 2 = 0, t c - 80. The orbit is plotted 
in configuration space x vs. t (left panel) and in phase space x 
vs. x (center panel), along with its corresponding power spec- 
trum (right panel) wherein the complexity -Ko.9 is provided 
as a measure of orbital chaoticity (see Sec. 1111 Ki . The four 
rows correspond to different noise strengths: (a) (|<5w|) = 0, 
(b) (|«5a;|) =0.001, (c) (|<ta|)=0.01, and (d) (|£w|)=0.1. 
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FIG. 2: Halo amplitude R H vs. t with T 2 =0, t c - 80, and 
(a) Ti = 0.05, (b) Ti = 0.1. The number of test particles is 
N = 10 6 . The four curves correspond to four different noise 
amplitudes. Blue curve with crosses: (|<5oj|) — 0. Red curve 
with asterisks: (|<5w|) = 0.001. Black curve with diamonds: 
(|5o>|) =0.01. Green curve with triangles: (|5w|)=0.1. 
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FIG. 3: Maximum halo amplitude max(Rii) reached over 
a duration t = 512 DE units vs. the logarithm of the test- 
particle population N with T2 = and t c = 80, and with (a) 
Ti =0.05, and (b) Ti =0.1. The four curves correspond to four 
different noise amplitudes. Blue curve with crosses: (|<5u>|) = 
0. Red curve with asterisks: (|5w|) = 0.001. Black curve 
with diamonds: (|5a>|) = 0.01. Green curve with triangles: 

<IH)=o.i. 
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FIG. 4: Maximum halo amplitude max(RH) vs. the loga- 
rithm of the autocorrelation time t c computed for TV = 10 4 test 
particles with Ti =0.05 and F2 = 0. Blue curve with crosses: 
(\Su>\) =0.002. Red curve with asterisks: (\6u>\) =0.01. Black 
curve with diamonds: (|&j|) =0.03. Green curve with trian- 
gles: (|5w|)=0.1. 
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FIG. 5: Percentage P{r>R) vs. R of test particles ( JV = 10 6 ) 
lying outside radius R at the end of the simulation (t = 512) for 
various noise strengths with (a) Ti =0.05, and (b) Fi =0.1. 
Fixed parameters are t c = 80 and F2 = 0. Blue curve with 
crosses: {\5uj\)—0. Red curve with asterisks: (|<5w|) = 0.001. 
Black curve with diamonds: (|<5w|) =0.01. Green curve with 
triangles: {\5ui\)=0.1. 
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FIG. 6: Plots of P(r > R) vs. i? in the presence of both 
the n — 1 and n = 2 modes (iV = 10 6 ) at the end of the 
simulation (t — 512) for fixed t c = 80, various values of mode 
amplitudes IY2, with noise strength (a) (|<5^|) =0.001, and 
(b) (\6u)\)=0.6l. Blue curve with crosses: I\ = 0.05, T 2 = 0.1. 
Red curve with asterisks: Ti = 0.1, T2 = 0.05. Black curve 
with diamonds: Fi =0.05, T2 ~0- Green curve with triangles: 

ri=o.i, r 2 =o. 
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FIG. 7: Evolution of four collections of 1600 test-particle 
orbits integrated over £ = 512 DE units (~60 to)- The collec- 
tions start (with zero initial particle velocity) at x = 0.31 ± 
0.002, x = 0.70 ± 0.002, x = 1.1 ± 0.002, and x = 1.41 ± 0.002 
with t c = 80, T 2 = 0, and: (a) I\ = 0, (\5u\) = 0, (b) I\ = 0, 
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FIG. 8: Halo distributions (N— 10 6 ) for radial orbits (dotted 
curves) and initially circular orbits (solid curves), (a) Halo 
amplitude R H (t) vs. t with ri=0.05, r 2 =0, and t c = 80. (b) 
Percentage P(r > R) of test particles lying outside radius R 
at the end of the simulation (t = 512). Blue curves: {\5uj\)=0. 
Red curves: (|<5u|) =0.001. Black curves: (|5o»|) = 0.01. Green 
curves: (|<5o>|)=0.1. 
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FIG. 9: Poincare sections for a set of 18 representative initial 
conditions integrated over £ = 2048 DE units (~250 to) with 
various noise strengths, and with i c = 80, Ti = 0.05, and 
T 2 = 0: (a) (|*w|> = 0, (b) (|«H> = 10- 6 , (c) (|$w|) = l(r B , (d) 
(|5w|) = 10- 4 , (e) (|<H) = 5 x 10_4 > and ( f ) (|<M) = 1 _3 - 
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FIG. 10: Poincare sections for a set of 18 representative initial 
conditions in an envelope-mismatched beam integrated over 
t = 2048 DE units (~250 to) with various noise strengths, 
and with M = 1.1118: (a) (\8lo\) = 0, (b) {\Su\) = 10~ 3 , and 
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